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Abstract 

Ph 

Various approaches aim to describe the recent analysis by the OPERA experiment, which 
indicates that neutrinos travel faster than the speed of light. We demonstrate that any 
such theoretical or experimental explanation must not destroy the complicated (nonlinear) 
structure of the proton waveform recovered in the neutrino signal. As one example, consider 
that only a fraction of the neutrinos travel faster than the speed of light, such as sterile 
neutrinos. We fit the OPERA data including this fraction as a free variable, assuming that 

t-h the OPERA result is correct. In our analysis, the best-fit values are 50% of the neutrinos 
i being superluminal and (v — c)/c = 4.5 • 10 -5 , where the neutrino velocity increases as 
£> the fraction of superluminal neutrinos decreases. The minimal fraction of superluminal 

^ neutrinos is found to be 17% (3a), which is constrained by the non-linearity of the proton 
c5 waveform. This minimal fraction challenges the hypothesis that only sterile neutrinos travel 
faster than the speed of light. In addition, we demonstrate that an experimental effect 
introducing a smearing between the proton waveform and neutrino signal, as expected for 
some systematical errors, is also limited by the shape of the waveform. Finally, we illustrate 
that even stronger constraints may be obtained from the recent analysis with a short-bunch 
beam, in spite of the low statistics. 
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1 Introduction 



Recently, the OPERA collaboration has measured the muon neutrino velocity, and (prelim- 
inarily) reported a difference between neutrino velocity and the speed of light (y — c)/c = 
2.37 • 1CT 5 |1|, which is compatible with an earlier MINOS result [2]. As far as the in- 
terpretation of this result is concerned, see Refs. [3]48] for general comments, Refs. (9 16 



for theoretical ideas indicating Lorentz invariance violation and Ref. 17 for earlier con- 
straints, Refs. 18-22 for ideas involving environment-dependent effects or new couplings, 
Refs. 23 30 for cosmological and astrophysical implications, and Ref. 31 32 for the dis- 
cussion of the analysis itself. In particular, any theory describing the superluminal motion 
has to face constraints such the bremsstrahlung issue discussed in Ref. [8). We do not give 
any theoretical explanation of the superluminal motion, in particular, we do not answer 
the question if neutrinos can travel super luminally or not. Instead, we discuss the addi- 
tional implications any theoretical or experimental description has to face, assuming that 
the result cannot be described by a trivial systematical error, such as an additional time 
correction not accounted for. More specifically, we show that the information contained in 
the non- linear shape of the proton waveform is important, where we give a theoretical and 
an experimental example. 

As one of the major observations, the OPERA results are obviously inconsistent with very 
stringent limits from SN1987A 33 -35 . This implies either a relatively strong energy depen- 



dence of the effect causing the superluminal motion of the neutrinos, see, e.g., Refs. 36 37 



or a flavor- or mass-eigenstate dependence, see Refs. [5j[8j[T2j[T6j|38j[39] for related discus- 
sions. A particular mass-eigenstate dependent case are sterile neutrinos, such as taking 
a shortcut through an extra dimension 40 44 , see also Ref. 45 for a similar approach. 



Especially, sterile neutrinos may be an elegant way to circumvent the argument in Ref. [8] , 
since they do not couple to the electroweak gauge bosons. On the other hand, the observa- 
tion of neutrino oscillations is, in general, incompatible with the hypothesis that only some 
of the active states travel faster than the speed of light 46 47 . However, this argument 



does not necessarily apply if only the sterile neutrinos are affected, since the incoherent 
limit, where active-sterile oscillations are averaged out, is perfectly consistent with short- 
baseline oscillation fits. In addition, note that neutrino oscillations have not been observed 
in OPERA yet [y T appearance can be described by flavor mixing), which means that even 
the propagation of the active neutrinos may be incoherent. In this letter, we use the most 
general assumption that only a fraction X of the observed neutrinos are superluminal. This 
fraction can be related to one or more of the active or sterile mass eigenstates, as it may be 
plausible for a tachyonic neutrino. In addition, we test the hypothesis if the signal may due 
to some additional smearing not accounted for in the analysis, which may be coming from 
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some unknown systematics. Finally, we comment on the recent test with a short-bunch 
wide-spacing beam. 



2 Neutrino propagation 



The OPERA experiment is sensitive to i/„ -> where the fraction of the mass eigenstate z/j 



in is |?7 M j| 2 . For the best-fit values from a recent global fit 48 , |t/ M i| 2 ~ 0.14, lU^] 2 — 0.35, 
and |f//i3| 2 — 0.51 if only active neutrinos exist, which means that mostly z/3 are produced in 
the OPERA beam. Since only are detected, the transition probability contains another 
factor of It/^l 2 . For instance, if the propagation of the neutrinos were incoherent, one would 



have 
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in the presence of active and sterile states. For the current best-fit values, — 0.02, 

It/^l 4 — 0.12, and |£7 M 3| 4 — 0.26, which means that v% clearly dominates. 

We assume that the fraction X of the muon neutrinos detected in OPERA is superluminal, 
traveling with a common speed (v — c)/c, and the fraction 1 — X is subluminal, traveling 
roughly with light speed. This assumption is very general in the sense that it does not 
depend on coherence (and, therefore, neutrino oscillations), which may only apply to a 
subset of the states - such as the active ones. For the special case of incoherent propagation, 
the fraction X of superluminal neutrinos can be written as (c/., Eq. (TO) 



E \u, 



/UJ I 



i, superluminal , . 

i, subluminal i, superluminal 



where, in general, all active and sterile states are summed over. If there are only active 
states, X ~ 0.05 for (only) u% superluminal, X ~ 0.30 for (only) v 2 superluminal, and 
X ~ 0.65 for (only) u 3 superluminal. If there are sterile states, for instance, a current 
short-baseline analysis of neutrinos oscillations 49 favors, in a 3+2 model, [U^l = 0.165 



and {Upsl = 0.148 at the best-fit (and also close to the incoherent limit, which may be 
applicable here), which means that |/7 M | 4 + l^sl 4 ^ 0.0012. Thus, if only the sterile states 
are superluminal, one has X ~ 0.003 in this simple model. In general, however, larger values 
in the percent range are potentially not excluded, whereas values X ^> 0.1 are inconsistent 



with other long-baseline disappearance data, see, e.g., Ref. 50 . 
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Figure 1 : Comparison of the measured neutrino interaction time distribution (data points) and the best- 
fit of the proton waveform (shifted by 1043.4ns; solid curve); figure reproduced from preliminary result in 
Ref. [IJ (v2). The dotted curve shows the proton waveform with a hypothetical (Gaussian) filter with a 
width of 200 ns applied. 

3 Method 



The measurement of (v — c)/c is based on the time delay 5t of the actual arrival time of the 
neutrinos compared to the expected arrival time assuming speed of light propagation, where 
a number of known corrections are taken into account and subtracted in the original analysis. 
For the derivation of this time delay, the shape of the proton waveform, measured by a beam 
current transformer, is matched to the shape of the neutrino signal. We have first reproduced 
the result (St = 57.8 ns) from Ref. [l], where we have only considered the statistical error 
and their first extraction. We show in Fig. [T] the corresponding comparison of the measured 
neutrino interaction time distribution (data points) and the proton waveform function (solid 
curve, shifted by 1043.4ns, which includes a known 985.6 ns correction), reproduced from 
Ref. jl] within the error from the resolution of the figures. Our analysis uses a Gaussian 
X 2 with a free normalization of the proton waveform to be fitted as another parameter. 
However, 5t and the proton waveform normalization have turned out to be uncorrelated, 
which means that the normalization has been determined and fixed for the further analysis. 
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One can easily see from Fig. [T] that not only the rising and falling edges between proton 
waveform and neutrino signal match very well, but also the nonlinear structure at the 
plateau. Thus, any effect modifying this structure will lead to a reduced goodness of fit. 

Our minimal \ 2 per d.o.f. is slightly larger (1.3) than the one of the OPERA analysis (1.1). 
We find that St = can be excluded at 5a, roughly consistent with the OPERA result, and 
our statistical error is slightly larger (about 11 ns, compared to 8 ns in the original analysis). 
These differences are expected from our rough approximation of the proton waveform, and 
from the different statistical method used here. We have also tested the impact of using 
more bins for the rising and falling edges of the the measured neutrino interaction time, 
without qualitative change. Since we prefer a self-consistent picture, we base our analysis on 
the waveform in Fig. [TJ Note that we do not consider systematical uncertainties quantified 
by the OPERA collaboration here, since they only add to the statistical uncertainty of the 
time of flight measurement in a trivial way, and cannot affect the shape of the waveform to 
the degree we discuss - as we will illustrate below. 

4 Fraction of superluminal neutrinos 

We assume that a fraction X of the detected neutrinos travels superluminally, and a fraction 
1 — X subluminally. We have first of all observed that a smaller X can be easily compensated 
for by a larger St. This near-degeneracy comes from the approximate linearity of the proton 
waveform at the rising and falling edges (c/., Fig. [TJ, which dominate the measurement. In 
the linear regime at the edges, one can show that the degeneracy 

X x5t~ 57.8 ns (3) 

holds. However, a too small value of X will lead to a change of the spectral shape in the 
non-linear regime, such as on top of the proton waveform, which limits the fraction X from 
below. 

We show our main result in Fig. [2] as a combined fit of St and X from the OPERA result 
(2 d.o.f.). The OPERA result corresponds to X = 1, with the best-fit marked by the red 
(upper) dot. It corresponds to the case of all states superluminal. If A < 1, one can 
clearly see the degeneracy with St, i.e., larger St can compensate for smaller X (cf., dotted 
curve). In fact, the 5a contour extends to very large values of St. This means that if X 
is marginalized over, St can be hardly constrained from the above. However, it can be 
clearly seen that St = can be excluded for any X at about 5cr, including the special case 
St = X = 0. The best-fit is found at St = 110 ns and X = 0.5, where the x 2 is relatively flat 
for X > 0.4. Therefore, only a lower bound for X can be determined, which is X > 0.17 at 
3a (1 d.o.f., since projected on vertical axis). 
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Figure 2: Combined fit of St (or (v — c)/c, on upper horizontal axis) and X from the OPERA result 
(2 d.o.f.). Here only the statistical error is included. The fractions X expected from Eq. ^ for specific 
superluminal mass eigenstates (All states, only, vi only) in the incoherent limit are also shown. The 
(preliminary) OPERA result corresponds to X = 1, with the best-fit marked by the red (upper) dot. The 
two-parameter best-fit is marked by the black (lower) dot. An additional constraint, which may be obtained 
from the test with a short-bunch wide-spacing beam, is shown by the dashed regions. The dotted curve 
represents the parameter degeneracy in Eq. ([3]). 

As far as the interpretation of our result is concerned, let us focus on Eq. ^ first, which 
assumes the incoherent limit. It can be read off in this simplistic approach that the OPERA 
result is consistent with only V2 or U3 being superluminal (see arrow markers in figure), 
whereas only v\ can be excluded. Of course, also combinations of the states or even all 
three states could be superluminal. If only sterile neutrinos are superluminal, such as by 
shortcuts though an extra dimension, the corresponding fraction X has to be X > 0.17 
in order to be compatible with the results at 3a. Note that in this case, the active states 
may still oscillate. Such a large contribution X of sterile neutrinos is, however, implausible 
with respect to short baseline data for sterile neutrinos in the electronvolt range. Therefore, 
theories with superluminal sterile neutrinos can probably be ruled out - unless X in the 
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detector is especially enhanced for the OPERA experiment parameters, and Eq. ^ is 
severely modified; see, e.g., Ref. 51 for a specific model, which demonstrates how to use 
this information in specific models. 

In the most recent OPERA paper [l] (v2), a test with a short-bunch wide-spacing beam was 
performed in order to trace back the individual neutrinos to specific proton bunches. In this 
test, no neutrinos with St = have been observed. This imposes an even stronger constraint 
on X, in spite of the low statistics of only 20 events. If, for instance, the distribution of 
St values is divided into two bins (split at St = 31 ns), at least 80% of the neutrinos have 
to be superluminal at 3a from a simple Possonian x 2 -test (including mis-identification of 
events, which can be estimated from the rms). In addition, St is directly constraint from 
the distribution of <5t-values, where the error can be obtained from the measured rms. In 
Fig. U these new constraints are shown as dashed regions for 3a and 6a. However, these 
bounds based on the assumption of a Gaussian distribution of the shape of the signal (St 
values), which is hard to confirm for the given low statistics. More general, it cannot be 
uniquely established that there is only one velocity component. 



5 Gaussian averaging of proton waveform 



In Ref. 32 , it was speculated that an additional smearing effect, not accounted for in the 
original analysis, may cause a deformation of the rising and falling edges of the neutrino 
waveform compared to the proton waveform, or even be implied in the proton waveform 
if not too large. Likewise, many systematical errors may not only cause a systematical 
time shift, such as the ones discussed by the OPERA collaboration, but also introduce 
some smearing. This includes also theoretical effects, such as non- negligible extensions of 
the neutrino wavepackets, see, e.g., Ref. [52j (Fig. 4). On the other hand, the shape of 
the proton waveform is convincincly well recovered in the neutrino signal, cf., Fig. [TJ even 
in the nonlinear regime. Therefore, the impact of such systematical errors cannot be too 
large. In order to quantify this statement for this particular class of systematical errors, we 
assume that such a smearing could be Gaussian in the simplest case, around a mean offset 
in St. For example, the dotted curve in Fig. [T] shows a smeared proton waveform where 
Gaussian smearing with a width of 200 ns has been used. Obviously, this smearing deforms 
the rising and falling edges. However, it also deforms the spectral features on top of the 
proton waveform, which will eventually reduce the goodness of fit. We leave the width a 
of such a Gaussian smearing as a free variable and fit it together with St in Fig. [3] As it 
can be read off from the figure, the best-fit St hardly changes in the presence of a possible 
smearing. Around the value a ~ 50 ns, where the rising and falling edges are deformed, 
the confidence level for the exclusion of St = slightly decreases. However, from the top of 
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Figure 3: Combined fit of St (or (v — c)/c, on upper horizontal axis) and a filter width a from the 
OPERA result (2 d.o.f.). Here only the statistical error is included. The (preliminary) OPERA proton 
waveform corresponds to a — > 0, with the best-fit marked by the red (lower) dot. The two-parameter 
best-fit is marked by the black (upper) dot. An additional constraint, which may be obtained from the test 
with a short-bunch wide-spacing beam, is shown by the dashed regions. 

the contours, o < 140 ns at the 3a CL (1 d.o.f.), where a significant contribution to the \ 2 
comes from the nonlinear part of the proton waveform. This implies that even asymmetric 
effects must be limited if they introduce an smearing over such timescales. In addition, 
from Fig. [3j 8t = can be excluded at the 4<r confidence level, even if a is marginalized 
over, which means that the OPERA result cannot have been compromised by such a large 
smearing effect (if Gaussian). In the meanwhile, as expected, no such effect has been found 
by the test with the short-bunch beam. In fact, if the shape of the distribution of events is 
assumed to be Gaussian, an additional constraint can be derived, shown as dashed curves in 
Fig. |3j In this case, the width of the smearing a is stronger limited. However, given the low 
statistics, the Gaussian nature of the distribution cannot yet be unambiguously established. 

Note that the class of systematics discussed here may also impact the measurement of X 
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discussed above, if present, since the measurement of X relies on the shape of the proton 
waveform. On the other hand, one can learn from Fig. [3] that smearing effects at the magni- 
tude of the systematical errors discussed by the OPERA collaboration (up to 8.3 ns) cannot 
distort the shape of the proton waveform significantly (in addition to introducing a sys- 
tematical shift). This may be different for the short-bunch measurement, where additional 
systematical errors dominate the measurement. 

6 Summary and conclusions 

We have demonstrated that theoretical or experimental interpretations of the OPERA result 
changing the shape between proton waveform and neutrino signal lead to a reduced goodness 
of fit. We have used two examples: the fraction of neutrinos being superluminal, and possible 
smearing effects between between proton waveform and neutrino signal. 

Assuming that only a fraction X of the neutrinos is superluminal, the best-fit is found to 
be at X = 0.5 with 5t = 110 ns, corresponding to (v — c)/c ~ 4.5- 10 -5 . This result is 
in fact closer to the MINOS result than the original OPERA result. However, fractions 
X > 0.35 are equally plausible. Most importantly, X can be limited from below X > 0.17 
(3<r), which means that superluminal sterile neutrinos will be challenged. On the other 
hand, only or z/3 superluminal may be consistent with the OPERA data from the time of 
flight measurement. This result will be, however, not compatible with neutrino oscillations 
if neutrino oscillations can be established in OPERA. Even stronger constraints can be 
obtained from the test with a short-bunch beam in spite of the low statistics, where no 
neutrinos with the speed of light have been seen. Therefore, in conclusion, the most plausible 
scenario to describe the OPERA result, if correct, is that all eigenstates travel faster than 
the speed of light, whereas flavor- or mass-eigenstate approaches to reconcile the OPERA 
and MINOS with SN1987A and other data are not promising. 

We have also demonstrated that there cannot be a significant smearing effect in the ex- 
trapolation from the proton to the neutrino waveform, as it may be expected for a class 
of systematical errors not accounted for. This result is supported by the short-bunch run, 
although the shape of the signal cannot yet be uniquely determined in this case. 
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